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ABSTRACT 



< 

O . We present evidence for a Galactic North-South asymmetry in the number 

^ ■ density and bulk velocity of solar neighborhood stars. The number density pro- 

file, which is derived from main-sequence stars in the Sloan Digital Sky Survey, 
shows a (North— South)/(North+South) deficit at \z\ ~ 400 pc and an excess at 
\z\ ~ 800 pc. The bulk velocity profile, which is derived from the Sloan Exten- 
^ . sion for Galactic Understanding and Exploration, shows a gradual trend across 

OO . the Galactic midplane as well as smaller-scale features. We speculate that the 

• . North-South asymmetry, which has the appearance of a wavelike perturbation, 

O ' is intrinsic to the disk. We explore the physics of this phenomenon through an 



analysis of the linearized Boltzmann and Poisson equations and through one- 
dimensional simulations. The perturbation may be excited by the passage of a 
satellite galaxy or dark matter subhalo through the Galactic disk, in which case 



;_i . we are witnessing a recent disk-heating event, 

c^ ■ 

Subject headings: Galaxy: kinematics and dynamics — solar neighborhood 

Disk galaxies are dynamic systems that can develop bars, spiral structure, and warps. 
They tidally disrupt satellite galaxies and dark matter subhalos, a process that leaves behind 
streams of stars and dark matter. Likewise, satellites continually perturb the disk of their 
host galaxy. 

Despite the existence of these time-dependent phenomena, the assumption that galaxies 
are in equilibrium (i.e., stationary in the potential) has been used extensively to interpret 
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certain types of observations. Best known, perhaps, are attempts to infer the surface den- 
sity and ver tical force near the Sun from the kinematics of stars in the solar neighborhood. 
OortI (I1932I ) pioneered this program and elements of his original method are present in vir- 
tually all studies of this type. Approaches to the so-called Oort problem generally involve 
an estimate of the distribution function for solar neighborhood stars, f {z, Vz), where z is 
the position relative to the midplane of the Galaxy and Vz is the vertical velocity. For 
stars in equilibrium, and not too far above the midplane, this distribution is a function of 
the effective vertical energy per unit mass, Ez = v'^/2 + "^{z), where \l/(z) is the gravita - 
tional potential in the solar ne ighborhood. Nunierou s rese archers including iBahcalll (1984 ). 
Kuiiken fc Gilmorel fll989al lb[). IPlvnn fc Fuch^ fll994h . and Iflolmberg fc Flvnnl J2000h have 
extended and improved upon Oort's method. 

In this Letter, we estimate the number density of stars with both photometric and 
astrome tric observations fr om SDSS-DR8, the Eighth Data Release of the Sloan Digital Sky 
Survey (JAihara et al.ll201ll ). We also use a smaller spectrosc opic data set from SEGUE, the 
Sloan Extension for Galactic Understanding and Exploration ( lYanny et al.ll2009l ). to estimate 
the bulk velocity and velocity dispersion. While our motivating interest was in revisiting 
the Oort problem, our aim in this letter is to demonstrate the existence of a North-South 
asymmetry in both the spatial density and velocity distribution of solar neighborhood stars. 
We explore the hypothesis that the asymmetry reflects a coherent wavelike perturbation, 
which is intrinsic to the disk. A theoretical discussion, together with a simple one- dimensional 
simulation, reveals that perturbations of this type are entirely natural. 

We begin with a discussion of the photometric data set. Since our primary interest is in 
differences in the stellar distribution North and South of the Galactic midplane, we choose 
matched regions from the survey above and below the Sun. Specifically, we select stars with 
Galactic coordinates 100° < / < 160° and 54° < |6| < 68° that reside within a perpendicular 
distance of 1 kpc from the line connecting the North and South Galactic poles. We determine 
dereddened apparent magnitudes in g, r, and i using the extinction maps of ISchlegel et al. 



(119981 ). In doing so , we assume that all of the stars are beyond the dust, which is thought to 
lie at 1^1 < 125 pc (JMarshall et al.ll2006l ). Only stars with 14 < r < 21 are included in order 
to avoid problems with saturation and photometric errors. We calculate absolute r-band 



magni tudes, M^, using r — i colors and the photometric parallax relation from iJuric et al. 



(120081 ) as well as g — i colors and the corresponding relation from llvezic et al.l (J2008al ) . The 
selection function, which depends on the geometry of the survey and our brightness cuts, 
is used to convert the observed number of stars to stellar number density with height. A 
detailed discussion of the selection function and various systematic effects will be discussed 
in a forthcoming publication. 
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In Figure 1, we plot the number density n(zobs), where Zohs = z — Zq and Zq is the 
vertical position of the Sun. We consider stars with 0.6 < r — i < 1.1 and compute distances 
using the photometric parallax relation from iJuric et al.l ( 120081 ). In this r — i range, we 



observe all stars (approximately 300K) with 0.2 kpc < l^obsl ^ 1.6 kpc. T he density profile 



show s the well-documented excess over a single exponential for \z\ > Ikpc (JGilmore fc Reid 
19831 : iReid fc Majewskilll993l : IJuric et al.l 120081 ). We fit a two- component model 



n{zohs) 



riQ I sech^ 



^obs + Zq 

2Hi 



+ /sech^ 



^obs + Zq 

2Ho 



T. 



using x^-minimization (jPress et al.l Il992l ) and find uq = (4.22 ± 0.04) x 10® kpc ^, Hi = 
(221 ± 3) pc, H2 = (582 ± 12) pc, / = 0.17 ± 0.01, and Zq = (15.4 ± 0.8) pc. The residuals 
A = (data — model) /model are largely an odd function in z. In the North, the data lie 
below the model for z ~ 400 pc and above the model for z ~ 800 pc while the situation is 
reversed in the South. We note that while equation 1 captures the general features of the 
measured n(zobs), the fit is relatively poor with x^ — 444 for 47 degrees of freedom. We 
strongly rule out a symmetric disk model, even when we assume some correlation of errors 
in adjacent color bins that might arise from the conversion of colors to distances via the 
photometric parallax relation. 

A consequence of the z-odd residuals is that the model parameters are sensitive to 
the range in z considered. In particular, if we fit only data within 500 pc, we find Zq ^ 
39 pc. Thus , our two determinations of Zq bracket the typical values found in the literature 
Jjuric et al.lbood l 



To explore the density structure further, we consider the North-South asymmetry, de- 
fined as A{z > 0) = {n{z) — n{—z)) / {n{z) + n{—z)). In Figure 2a, we show A for a series 
of color bins from r — i = 0.4 to r — z = 1.4. In constructing this figure, we use Zq = 15 pc. 
Note that the probed range in z depends on color; regions close to the midplane are probed 
by the redder bins and vice versa. Evidently, A shows little color dependence. In Figure 2b, 
we show A for the r — i window used in Figure 1. Also shown is the asymmetry obtained 
if distances for these stars are computed using the g — i photometric parallax relation. In 
Figure 2c, we compute the composite asymmetry over all color bands in Figure 2a as well 
as the composite asymmetry found for stars with 1.6 < g — i < 2.6. We also show the 
asymmetry when we assume Zq = 39 pc. 

The statistical significance of the North-South asymmetry may be tested using the 
Kolmogorov- Smirnoff (K-S) test where the statistic D is defined as the maximum difference 
between fractional cumulative distributions of stars North and South of the Galactic mid- 
plane. We carried out a K-S test for subsamples with different numbers of stars A^ and 
found that D ~ 0.035 - 0.055 for A^ > 9000. The probability that the stars North and South 
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of th e midplane are drawn from the same distribution is P ~ exp {—D'^N) fjPress et al. 
19921 ). Thus, with only a few percent of the full data set we rule out the null hypothesis 
(North-South symmetry distribution) at the 10~^ level. For the full data set, P is vanishingly 
small. 

In Figure 3 we show the two-dimensional asymmetry, A2D (R, z ), where R is the G alac- 
tocentric radius. This figure may be compared with Figure 26 from lJuric et al.l (120081 ). who 
also studied SDSS stars. The asymmetry described in this letter appears to be evident in 
the bottom right panel of their figure. 

We now consider the spectroscopic sample, which comprises some 11,000 SEGUE stars. 
The (/, 6)-footprint for this sample is slightly larger than the footprint for the photomet- 
ric sample, though the (/,6) coverage is not uniform (JYanny et al.l 120091 ). We inclu d e star s 
with 1.0 < f7 — 2 < 1.9 that satisfy the stellar surface gravity cut from iLee et al.l (120081 ) . 
which restricts our sample to dwarfs. Proper motions and radial veloc ities are converted to a 



v^u 



yi 



Vz) coordinate system using the standard transformation laws (jjohnson fc Soderblom 



19871 ). Errors in v^ are typically 6 — llkms"^ and result from uncertainties in proper mo- 
tions, radial velocities and distances. We consider stars with \vz\ < 125kms~^ so as to avoid 
contamination from high- velocity outliers. In Figure 4a we see a trend in (vz) across the 
midplane as well as small scale undulations. These results suggest that there is a coherent 
motion of solar neighbo rhood stars away from the Galactic midplane. We note that recently 
Bochanski et al.l (I2OIII ) found evidence for variations in (vz) with M-star color, which is cor- 
related with z. Figure 4c shows an increase in cr^ with \z\, which may be attributed to an 
increase in the fraction of kinematically "hotter" thick disk stars as one moves away from 
the midplane. 



Various systematic effects can produce trends in (vz). iBond et al.l (120 lOl ) analysed SDSS- 
DR7 data and reported a difference of ~ 8 kms~^ between the bulk motions North and South 
of the midplane; they attributed this difference to a systematic radial velocity error for M 
dwarfs possibly due to the limited set of velocity template spectra used in the SDSS cross 
correlation. To explore this issue, we consider red and blue subsamples. The results, shown in 
Figure 4b, do indeed indicate a stronger trend in the red subsample, though a trend and small 
scale features are still evident in the blue subsample. Another example is a systematic offset 
in distanc es, which would enter (f^) through the proper motions, especially at mid-Galactic 
latitudes (ISchonrich et al.ll2012l ). though a large offset (~ 30%) is required to explain the 
trend seen in Figure 4a in toto. Clearly, a definitive statement on the z-dependence of (f^) 
will require improvements in radial velocity and distance measurements along with better 
coverage toward the Galactic poles. 



A primary goal of the SDSS and SEGUE projects is to discover new structure in the 
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Galaxy's stellar halo. iJuric et al.l ( 120081 ). for example, detected several new overdensities in 
the thick disk and stellar halo. These features are usually interpreted as stellar debris from a 
tidally disrupted satellite. While it is tempting to seek a similar explanation for the features 
described in this letter we explore a different interpretation: the structures seen in both the 
photometric and spectroscopic data represent coherent wavelike perturbations in the Galactic 
disk. 

Evidence for this interpretation can be found in the preceding figures. The model 
residuals (Figure 1) and North-South asymmetry (Figure 2) extend above and below the 
Galactic midplane and exhibit a wavelike structure. Furthermore, the structures appear to 
run parallel to the Galactic midplane, as seen in Figure 3. In addition, the bulk velocity 
exhibits a trend running through the midplane and smaller-scale fluctuations. Finally, and 
perhaps most importantly, the asymmetry appears to be largely independent of color (Figure 
2a). A strong color dependence would have suggested that the as ymmetry was due to 



different stellar populations North and South of the midplane. Indeed, ISchlafly et al.l (120101 ) 
note that F turn-off stars in the South appear to be systematically redder in (7 — r by 0.02 
mags than in the North. However, these stars are typically located 5 — lOkpc from the 
Galactic plane, and therefore well beyond the extent of the thin or thick disk stars studied 
here. 

We now demonstrate that these features naturally arise in a perturbed self-gravitating 
stellar disk. For simplicity, we ignore interstellar matter and model the Galactic disk in the 
solar neighborhood as a coUisionless system of stars. We further assume that the disk is 
self-gravitating, though it is easy to include additional components to the gravitational field. 
Finally, since we work close to the Galactic plane, we assume that we can decouple vertical 
motions of disk stars from radial and azimuthal motions. 

We consider small perturbations to the equilibrium solution for a single- component, 
self-gravitating, isothermal plane. The distribution function, density, and potential for the 
equilibrium solution are 

foiz,v.) = -^^^2^^V{-EJ^')^ (2) 



^o(^) = 2a' In (cosh ( — 11 , po{z) = posech' ( — ) , (3) 



where po = (j'^/SirGz'^ ( (ISpitzeiill942l ). see, also iBinney fc Tremaind (120081 ) . problem 4.21). 



The linearized Boltzmann and Poisson equations are 



^ + ^ ^ _ ^^ _ ^^ = (4) 

dt ^ dz dz dv, dz dv^ 
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and 

-Q^ = 4vrGpi , (5) 

where pi {z, t) = J dv^fi {z, v^, t). Equation [5] implies that the vertical force perturbation, 
Fi = —d'^i/dz, is proportional to the perturbed surface density, as defined from the Galactic 
midplane, i.e., Fi{z) = AnG (T,i{z) — Si(0)) where Si (z) = J^^dz pi{z). Moreover, we can 
integrate equation H] over Vz to obtain the linearized continuity equation 

where {pv)i = J dvzfiVz- The mean velocity is therefore related to the surface density and 
vertical force: 

{pv\ - ^ - ArrG^ , (7) 

which implies a connection between the vertical perturbations considered here and the Oort 
problem. 

Equations H] and [5] do not admit simple plane wave solutions since the unperturbed 
model depends on z. Nevertheless, the standard Jeans stability analysis suggests that waves 
are described by the dispersion relation 

*=~^ + ^, (8) 

where k and u are the wavenumber and angular frequency, respectively. A lower bound on 
the Jeans length, the minimum length scale for gravitationally unstable modes, is obtained 
by setting u = and p and a equal to values characteristic of the Galactic midplane. We 
find \j = 27r /kj = (na'^/Gp) ^ 2kpc. Since Xj is much larger than the characteristic 
scale height of the disk, we conclude that vertical perturbations are stable and behave like 
pressure-supported waves. Equation [8], together with the fact that p/a"^ decreases with \z\ 
implies that the wavelength for these perturbations will be somewhat less than Xj and will 
increase as one moves away from the midplane. These conclusions are consistent with the 
asymmetry shown in Figure 2. 

Since the equilibrium distribution function is even in both z and Vz, the characteristic 
modes of the perturbed system have definite parity in both phase space coordinates. Any 
perturbation can be modeled by displacing stars in position and velocity, that is, by replacing 
the unperturbed phase space positions, zq and Vq with z = zq + Zi{zq) and Vz = Vq + Vi{zo). 
For the simplest perturbation, all stars are displaced in the same direction. We refer to 
this type of perturbation as an m = mode since there are no nodes in the displacement 
functions Zi and vi. The m = mode amounts to a local displacement of the stellar disk 
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as a whole and affects any determination of Zq and Vz^q from solar neighborhood data. 
(For observational evidence of m = modes or midplane displacement in the molecular gas 
distribution of the Galaxy, see iNakanishi &: Sofud ( 120061 ) .) 



The next simplest perturbation is an m = 1 (breathing) mode where displacements and 
peculiar velocities have opposite signs above and below the midplane. However, the density 
perturbation is an even function in z and therefore n{z) alone cannot reveal the presence of 
an 771 = 1 perturbation. However, the m = 1 mode would show up as a North-South trend 
in the bulk velocity, which is consistent with what we see in Figure 4. 

The North-South asymmetry seen in Figure 2 is then likely a combination oi m = 
2, 4, . . . modes. The local extrema in A at 400 and 800 pc (and possibly 1200 pc) imply 
nodes in zi at the corresponding positions in z. These higher order modes may also be 
present in the velocity data. 

We have numerically solved for the evolution of a one-dimensional collisionless N-body 
system where each "particle" corresponds to an infinite plane, which moves under the influ- 
ence of the gravitational field generated by the other planes. The gravitational force on a 
given plane is proportional to the difference between the number of p lanes to the North and 



the nu mber to the South. (Similar simulations were performed by ISanchez-Salcedo et al. 



(120 111 ) in order to study vertical equilibrium solutions of galactic disks.) 



We chose an equilibrium model given by equations [2] and [3] with a = 31kms^^ and 
Zd = 500 pc. Thes e parameters imply a surfac e density at 1.1 kpc of 71 Mq pc~^, which is the 
value obtained by lKuijken fc Gilmord (jl989al ). The initial perturbation is chosen to roughly 
match the results found in Figures 1 and 4. The evolution of the system is shown in Figure 
5. Ill The simulations illustrate the interplay between density and velocity perturbations. 
Waves appear to reflect off the low-density regions at high z, a phenomenon reminiscent of 
waves propagating in the Sun, and decay over a period of 200 — 300 Myr. 

To conclude, we speculate on the origin of vertical perturbations in galactic disks. An 
isolated, self-gravitating disk can develop non-ax isymmetric structure , such as sp iral arms 
and b ars, through gravitational instabilities (see iBinney fc Tremaind (|2008[ ) and ISellwood 
(I2OIOI )). On the other hand, the vertical modes are stable and must be excited by some 
agent. Of course, stellar disks of even relatively isolated galaxies, such as the Milky Way, 
interact with satellite galaxies and, presumably, halo substructu re. There is an extensiv e 



body of literature on disk heating, dating back to the work of iLacev fc Ostriker 



who considered disk heating by supermassive black holes, and IVelazquez fc Whitd ( 1l999l ) 



1985 ) 



^See www.astro.queensu.ca/~widrow/OortExplorers.htnil for an animation of the simulation. 



who considered disk heating by satelhte galaxies. In addition, iGauthier et al.l (|2006[ ) and 
Dubinski et al.l ( 120081 ) showed that satelhtes and halo substructure cai i provoke spir a l-type 
structures in a stellar disk and trigger bar formation. More recently, iPurcell et al.l (1201 if ) 
proposed that the Sagittarius dwarf galaxy was responsible for setting up the particular 
pattern of spiral structure seen in the Milky Way. Our conjecture is that the vertical waves 
are another manifestation of disk heating and disk-substructure interactions. 

The observational, theoretical, and numerical analyses presented in this paper can be 
extended to incorporate variations in the disk plane and provide a more complete pic- 
ture of the perturbed disk. On the observat ional front, f uture experiments, such as GAIA 
( iPerryman et al.ll2001t IWilkinson et al.ll2005[ ) and LSST ( llvezic et al.ll2008bl ). which aim to 
measure distances and velocities for some one billion stars, will provide a picture of the phase 
space distribution of the Galaxy in unprecedented detail. These studies may well reveal a 
rich structure of waves in the disk, and thereby open a new window into the dynamics of the 
Galaxy. 
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Fig. 1. — Number density n as a function of distance from the Sun 2;obs- Black points are 
the data. The magenta curve is our model fit. The dotted blue curve is the contribution 
from the thin disk; the dashed magenta curve is the contribution from the thick disk. Lower 
panel shows residuals: A = (data— model) /model. 
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Fig. 2. — North-South asymmetry parameter A as a function of z. Top panel shows A{z) for 
0.2 mag color bins: black square - r — i = {1.2 — 1.4}; red filled triangle -r — z = {l — 1.2}; 
green open triangle - r — i = {0.8 — 1}; blue filled circle - r — i = {0.6 — 0.8}; magenta open 
circle - r — i = {0.4 — 0.6}. Middle panel gives the average A for the color bins and z-range 
used in Figure 1. B lack points were generated using the photometric parallax relation from 
Ivezic et al.l (J2008al ). Bottom panel gives average A over the full range in z and r — i. Red 
triangles assume Zq = 39 pc instead of Zq = 15 pc. 
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Fig. 3. — North-South asymmetry A2D as a function of z and the Galactocentric distance 
R. We assume 8 kpc for the Sun-Galactocenter distance. 
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Fig. 4. — Bulk velocity (vz) and velocity dispersion cr^ as a function of z in units of kms^^. 
Top panel shows the bulk velocity as a function of z for the entire spectroscopic sample. 
Middle panel shows the bulk velocity profile for the "red" subsample {g — r > 1) and "blue" 
subsample {g — r < 1). The peculiar motion of the Sun (f^,© = 7.2kms~^ ( iDehnen fc Binney 
19981 )) has been subtracted from (f^). Bottom panel shows the velocity dispersion. 
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Fig. 5. — Density, bulk velocity and vertical force as a function of z for the initial conditions 
and for three time frames from the one-dimensional simulation described in the text. Top 
panel shows A = {p{z) — po) /Po where po is the unperturbed density. Middle panel shows 
the bulk velocity. Bottom panel shows the difference between the vertical force for the 
perturbed and unperturbed disk. Blue points represent results from our analysis of the 
SDSS and SEGUE data. 



